Vedeți discuțiile, statisticile și profilurile de autori pentru această publicație la: Vizualizări continue relaxate ale sistemelor cu evenimente discrete: considerații privind diagramele Forrester și rețelele Petri Lucrarea Conferinței ■ Noiembrie DOI: /ICSMC - Sursa: IEEE Xplore CITAȚII CITEȘTE autori, printre care: Universidad de La Rioja (Spania) PUBLICAȚII CITAȚII Universitatea din Zaragoza PUBLICAȚII CITAȚII VEZI PROFIL VEZI PROFIL Universitatea din Zaragoza PUBLICAȚII CITAȚII VEZI PROFIL Unii dintre autorii acestei publicații lucrează și la aceste proiecte conexe: Modelare și simulare în industria alimentară pentru sprijin decizional Model de dinamică a sistemului pentru evaluarea și simularea amprentei de apă a lanțului de aprovizionare cu bioetanol Conținutul AU care urmează acestei pagini a fost încărcat de la iulie Utilizatorul a solicitat îmbunătățirea fișierului descărcat IEEE International Conference on Systems, Man and Cibernetica Vizualizări relaxate continue ale sistemelor de evenimente discrete: considerații asupra diagramelor Forrester și rețelelor Petri* E Jiménez'* , J Júlvez *' ), L Recalde** , M Silva'' ' (* ) Dept Ingeniería Eléctrica Área de Ingeniería de Sisteme și Automática Universitatea de La Rioja Luis de Ulloa , E- , Logroño (Spania) e-mail: emilio jimenezm@die unirioja es (* ) Departamentul de Informatică și Inginerie a Sistemelor Universitatea din Zaragoza Maria de Luna CP- Zaragoza e-mail: {silva,lrecalde,julvez}@unizar en Rezumat - Rețelele Petri (PN) constituie o paradigmă formală pentru sistemele discrete Unele modele discrete pot fi relaxate în modele continue Rețelele Petri semantice infinite de servere (ISSCPNs) este una dintre cele mai relevante interpretări cronometrate ale Continuons PN ISSCPN-urile pot fi văzute ca sisteme liniare pe bucăți Diagramele Forrester (FD) sunt instrumente de modelare specifice din cadrul System Dynamics, o metodologie pentru analiza sistemelor continue complexe Această lucrare explorează și compară puterea de modelare a ambelor formaliste, ÌSSSCPN și FD La baza acestei lucrări stau pregãtirile comparative axate pe formalizatii si pe sistemele pozitive si compartimentale Comparația este completată ținând cont de interpretarea sistemelor de ecuații diferențiale liniare ordinare (LODES), de întârzierile de informare și de unele considerații metodologice ISSCPN-urile permit modelarea oricăror CODURI atunci când există limite superioare și inferioare cunoscute ale variabilelor de stare Prin urmare pot fi modelate sisteme cu comportament ciclic sau întârzieri în informație Cuvinte cheie: rețele Petri continue, diagrame Forrester, relaxarea sistemelor dinamice de evenimente discrete, sisteme pozitive, putere expresivă Introducere PN-urile constituie o familie binecunoscută de formalism dinamic al evenimentelor discrete peste naturais nenegativ Deși modelele PN sunt inițial modele de evenimente discrete, relaxarea lor prin continuizare le transformă în modele continue Cu prețul pierderii anumitor posibilități de analiză, aceasta permite obținerea unor avantaje, cum ar fi evitarea problemei de explozie a stării inerente Sistemelor discrete și profitarea teoriei extinse despre Sistemele dinamice continue Deși nu toate Sistemele PN permit o continuare „rezonabilă” [ ], această relaxare este posibilă în multe cazuri practice, conducând la un formalism în timp continuu: continui PN Diferite interpretări cronometrate duc la ' - - - / /$ © IEEE о Sprijinit prin grant de la DGA refBl / prevăzută cu un marcaj inițial peste P, Шо- Pre și Post reprezintă structura statică a modelului, față de care poate fi dedusă matricea fluxului de simboluri, C = Post - Pre O structură a rețelei Petri poate fi reprezentată și ca un graf bipartit direcționat, în care locurile p sunt de obicei reprezentate ca cercuri și tranzițiile t ca bare Marcajul se modifică prin declanșarea tranzițiilor Pornind de la un Sistem PN, se poate scrie o ecuație de stare (sau fundamentală): m = mo + С · σ, unde σ e N și me Næ Setul de stări accesibile ale unui sistem PN discret poate deveni cu ușurință extrem de mare Aceasta este cunoscută sub numele de problema exploziei de stat O modalitate de a încerca să depășim această problemă este de a continua sistemul, ceea ce permite utilizarea diferitelor matematici! instrumente (geometrie convexă și programare liniară, ecuații diferențiale ) PN-urile pot fi interpretate cu marcarea în reali nenegative (PN-uri continue) pentru a modela aproximarea continuelor relaxate a modelului discret Tragerea se modifică în același mod, adică se activează o tranziție t la m iff V pe 't, m[p]> Gradul său de activare este definit ca enab[r] = min [m[p]/Pre[p,z]) Într-un PN continuu putem lua în considerare derivata ecuației stării în raport cu timpul Astfel obținem că m = C ■ σ, plus condiția inițială m( ) = mo Să cali f = σ, deoarece reprezintă fluxul prin tranziții Dacă f(t), unde τ reprezintă timpul, este definit printr-o extensie interpretativă, atunci evoluția temporizată a PN continuu poate fi obținută Au fost definite diferite semantici de ardere pentru rețele Petri continue [ ], [ ] Semantica serverelor infinite va fi folosită în această lucrare Sub această semantică de ardere curgerea printr-o tranziție r; este definită ca П^}=Х[г,-]-епаЬ[?Д unde X[r,J este o constantă reală pozitivă reprezentând viteza internă a tranziției r; Definiții diagrame Forrester SD constituie o metodologie de modelare care permite sistematizarea creării de modele continue bazate pe Sisteme de ecuații diferențiale neliniare multivariabile dependente de timp FD sunt instrumente specifice de modelare din cadrul SD [ ] care oferă o reprezentare grafică a Sistemelor dinamice (vezi Figura ), modelând cantitativ relațiile dintre părți prin intermediul unor simboluri, care corespund unei interpretări hidrodinamice a Sistemului Sursă și siak (cloiKh) LeveJ variabilă (Stoc) Debit variabil (supapă) Întârziere (material sau informații) Variabila auxiliara Variabila exogcnùui Parametru (constante) Material cbannel (sau flux) Canal de informare! Nivelurile corespund variabilelor de stare din teoria sistemelor Ele reprezintă variabilele a căror evoluție este mai semnificativă pentru studiul Sistemului Nivelurile acumulează „material” din canalele de material, care sunt controlate de supape (variabile de debit) Supapele definesc comportamentul sistemului, deoarece determină viteza fluxului de material (prin canalele de material) conform unui set de ecuații asociate Ecuațiile depind de informațiile pe care supapele le primesc de la Sistem (niveluri, variabile auxiliare și parametri) și de la mediu (variabile exogene) Informația este transmisă instantaneu prin canale de informare Variabilele auxiliare corespund pașilor intermediari în calculul funcțiilor asociate supapelor, dar pot fi oricând eliminate Norii reprezintă surse și chiuvete, adică o cantitate nedeterminată (infinită) de material, iar parametrii sunt valori constante ale Sistemului Interacțiunea sistemului cu exteriorul este reprezentată de variabilele exogene, care au o evoluție care se presupune că este independentă de evoluția sistemului Întârzierile pot afecta materialul sau transmiterea informațiilor, dar în ambele cazuri nu introduc mai multă capacitate de descriere, deoarece corespund doar unei notări compacte a elementelor care produc aceste întârzieri Un model bazat pe ecuații diferențiale este derivat direct dintr-un FD, iar interesul analogiei hidrodinamice este că indică faptul că un model FD este echivalent cu un ODES (eventual neliniar) și invers Rețele Petri și sisteme de ecuații diferențiale obișnuite În secțiunea anterioară s-a văzut că evoluția unui ISSCPN este descrisă de sistem: m(T)= C ■ f(î) f(T)[ti]=A[ti] enab(T)[ti] m( ) = lună Astfel, o rețea Petri continuă sub semantică infinită de servere devine un sistem liniar pieewise Comutarea între două Sisteme liniare este declanșată de o schimbare a locului care dă minimul în expresie pentru gradul de activare Despre pozitivitate În linii mari, sistemele pozitive sunt sisteme ale căror variabile de stare iau doar valori nenegative Un sistem pozitiv păstrează automat non-negativitatea variabilelor de stare, adică „dacă se adaugă constrângeri de non-negativitate asupra stării, acestea sunt redundante Figura : Elementele Diagramelor Forrester Mai formal, fie £ ( ) un sistem liniar: x(i)=Âx(t) + Hu(t) ( ) Definiția [ ] Se spune că J este pozitiv dacă pentru fiecare stare inițială nenegativă și pentru fiecare intrare nenegativă starea sa este nenegativă Atunci orthandul pozitiv SK* este o mulțime invariantă nenegativă Dacă B= , se spune că Sistemul este necontrolat sau neforțat Rețineți că pozitivitatea în sistemele liniare poate depinde atât de baza intrării, cât și de baza spațiului de stare Un sistem nepozitiv poate fi transformat într-un alt sistem pozitiv echivalent printr-o schimbare de bază în spațiul stărilor Acesta este motivul pentru care unii autori definesc Sistemele pozitive cerând existența unei mulțimi invariante (fără a cere, totuși, ca un astfel de set invariant să fie orthandul pozitiv) Teorema [ ] Un sistem liniar ( ) este pozitiv, dacă A este o matrice Metzler și B este nenegativă (o matrice/vector este nenegativă dacă toate elementele sale sunt nenegative și o matrice pătrată este Metzler dacă elementele nediagonale sunt nenegativ) Conform definiției , ISSCPN-urile sunt sisteme pozitive (faptul că fluxul unei tranziții este proporțional cu gradul său de activare asigură nonnegativitatea marcajului) Cu toate acestea, matricele A,· ale sistemelor liniare care guvernează evoluția rețelei (recuperând că un ISSCPN este un sistem liniar pieewise) pot fi matrici non-Metzler Într-un ISSCPN comutarea între sisteme liniare este declanșată de o schimbare a locului care dă minim în exprimarea gradului de activare Evoluția sistemului net din figura este condusă de un sistem liniar cu matricea A] dacă xi ci-S)/(c c )) ) ' Deci, dacă este reprezentat Depozitul versus angajare, deși vânzările sunt constante, apare un comportament ciclic, cu parametrii indicați în Figura Ca curiozitate, acest tip de comportament (ciclic chiar și cu intrări continue) a fost la originea studiilor lui Forrester , care au fost sursa dinamicii sistemului Figura Comportamentul sistemului descris în Figura Dacă St( )=DSt și E( )=S/cb atunci sistemul este stabil De asemenea, este important de subliniat faptul că acest sistem are sens fizic numai atunci când nivelurile (elementele stocate și numărul de angajați) au valori pozitive, dar sistemul de ecuații diferențiale nu este pozitiv, iar angajarea și stocarea negative pot fi atinse pentru unele condiții inițiale Prin urmare, constrângerea pentru nenegativitate (St, E> ) trebuie inclusă suplimentar pentru a obține un model corect Sistemul cu constrângere de non-negativitate poate fi modelat cu PN continue (Figura ), dar trebuie luat în considerare faptul că C, care poate fi pozitiv sau negativ, trebuie implementat ca o combinație între un flux de contracte noi și un flux de concedieri, ambele pozitive Rețineți că au fost utilizate două locuri (cele cu marcaj unitar) pentru a obține un flux constant cu ISS, iar arcuri de control au fost necesare pentru a selecta în mod explicit locurile care furnizează informațiile tranzițiilor cu sincronizări Sistemul va fi descris prin ( ) ori de câte ori St>l/kx și E>St/k„ Recaliem că kc reprezintă o constantă finită cât de mare ne dorim (în cele din urmă tinde spre infinit) Figura ISSCPN echivalent cu sistemul din Figura limitat la valori pozitive De asemenea, este important de menționat că o valoare adecvată a lui ka, în PN depinde de valorile minime pe care le pot atinge St și E (și apoi de marcajul inițial) De exemplu, dacă d=c =l, St( )= , DSt= și E( )=S= , atunci k„ nu trebuie să fie mai mare de Figura arată simultan evoluția Sistem constrâns (modelat cu ISSCPN sau FD cu constrângeri) și cel nerestricționat (modelat cu ODES sau FD fără constrângeri) Ambele sunt similare de la starea inițială până la prima intersecție cu axa orizontală (punctul a în grafic) Figura a prezintă St versus E, iar Figura b arată evoluția temporală a variabilelor de stare Evoluția angajaților vs stocare - Depozitare (Sf) ♦ці" ■ ■ - ——— Ѳ Timp Figura : Evoluția sistemului cu ci=cj=l, St( )=E( )= ; DSt= , C=V= din modelele constrânse (linie întreruptă) și fără constrângere (linie întreruptă) Întârzieri de informare în PN Modelele cu întârzieri de material și informații sunt utilizate frecvent în FD În ambele cazuri, întârzierile sunt modelate fie ca un sistem de prim ordin, fie ca mai multe în serie (de obicei trei), dar nu sunt utilizate elemente de modelare suplimentare În PN întârzierile materiale sunt adesea folosite, similar întârzierilor materiale în FD, dar întârzierile similare pot fi, de asemenea, modelate cu informații Tabelul prezintă ambele tipuri de întârzieri în FD și în ISSCPN-urile lor echivalente Rețineți că dinamica (ecuația diferențială) este similară în ambele întârzieri (un sistem de ordinul întâi), dar întârzierea materialului este aplicată unui flux (canale/arce de material), în timp ce întârzierea informației este aplicat la informaţia unei variabile de stare (canale de informare de la un nivel ! arcuri de control cu informaţia unui loc) Tabelul Schiță, ecuație diferențială, FD și ISSCPN-uri în întârzieri materiale și informaționale Să presupunem un alt sistem de stocare Vânzările (care sunt retururi atunci când sunt negative) sunt proporționale (cu constanta de proporționalitate k) la diferența dintre referința de stocare (r) și numărul de elemente stocate (y) Figura prezintă sistemul modelat cu un FD și, respectiv, cu un ISSCPN Figura a) DF care modelează un sistem de stocare simplu, b) PN echivalent al acestuia Acest sistem urmează ecuația diferențială dy/dt = ku = - k ■ y + k · r A cărui soluție este y = (y - r) ■ e'kl + r, unde y este stocarea inițială Rețineți că aceste ecuații sunt aceleași cu cele ale sistemului tipic de control al temperaturii (y), cu o buclă de feedback de referință de temperatură (r) și eroare (u) Să presupunem acum o întârziere în bucla de feedback informațional, adică informația de eroare u este întârziată cu /a unități de timp, iar x este variabila care reprezintă semnalul întârziat Dacă este folosită o aproximare de ordinul întâi pentru a modela întârzierea, ecuațiile diferențiale care conduc evoluția sistemului sunt x = -axay + a- r ÿ = к · x Și Figura prezintă modelele cu FD și PN-uri, respectiv Figura Modelele FD și ISSCPN ale unui sistem de stocare cu o întârziere de ordinul întâi în bucla de feedback informațional Rețineți că, din cauza pozitivității PN-urilor, PN din Figura nu este valabil pentru valorile negative ale variabilelor de stare (și apoi ale lui x) Și x devine negativ ori de câte ori y(O)>r Pentru a evita această problemă, se poate face o schimbare de referință a originii când x este mărginit inferioară: x = x'-lim unde lim” este o constantă care reprezintă o limită inferioară a lui x Ecuațiile rezultă atunci x'=-a x'-a-у+ar + a-tim ( ) ÿ = kx'-klim Figura arată PN care modelează sistemul cu modificarea originii referencé (sau ODES în ( )) Figura Modelul ISSCPN al Sistemului de stocare cu o modificare a variabilei de stare: x - x'-lim Cu toate acestea, pentru acest exemplu, există o soluție mai bună la problema de non-negativitate care poate fi aplicată chiar dacă x este nemărginit Eroarea и este o variabilă nepozitivă, dar corespunzătoare diferenței dintre două variabile pozitive: temperatura y și referința r Atunci x poate fi modelat prin două variabile, xr corespunzând reférencé întârziat și xy corespunzător temperaturii întârziate, așa cum se poate observa în Figura IL Acest model este corect pentru orice valoare inițială (nenegativă) și nu este necesar să se cunoașteți o limită inferioară a lui x Figura Modelarea PN a unei întârzieri de informare de prim ordin într-un sistem de stocare Cu această aproximare de ordinul întâi a unei întârzieri, matricea A a sistemului și valorile sale proprii sunt: -a -a к —±—Va"- -ak Aceste valori proprii au întotdeauna o parte negativă reală, iar apoi sistemul neforțat este stabil pentru orice timp de întârziere finit Este bine cunoscut faptul că o întârziere de aproximare de ordin al treilea este mai asemănătoare cu o întârziere pură și poate provoca comportamente instabile (poli cu parte reală pozitivă) Figura prezintă FD-ul sistemului de întârziere de ordin a treia și ODES-ul corespunzător dxi/dt = За(г-у-Хі) dxj/dt = a(X|—Xi) dxj/dt = a-(x -x,d dy /dt = xj Figura Modelarea DF a unui sistem de stocare cu întârziere de informare și ODES echivalent Acest ODES devine instabil când întârzierea crește, iar apoi apar stări negative Totuși, variabila de stare, care reprezintă elementele stocate, are o natură pozitivă Apoi sistemul trebuie constrâns la stări nenegative (y> ) pentru a obține un model corect Acest ODES (sau FD) cu constrângere nenegativă are același comportament ca modelul ISSCPN, care este prezentat în Figura Figura PN care modelează o întârziere a ordinului de sete Figura arată simultan, așa cum a făcut Figura în ultimul exemplu, evoluția sistemului constrâns (PN-uri și FD cu constrângeri explicite) și a celui nerestricționat (ODES și FD fără constrângeri explicite) Ambele comportamente sunt din nou similare de la starea inițială până la prima intersecție cu axa orizontală (punctul a în grafic) Cu PN-uri, prin schimbarea referinței, se poate simula o evoluție asemănătoare cu cea a sistemului nerestricționat de la punctul inițial până la intersecția cu orice axă paralelă cu axa orizontală Și această axă paralelă poate fi atât de negativă pe cât ne dorim Acest lucru este interesant din punctul de vedere al puterii expresive a ISSCPN-urilor Cu toate acestea, PN-urile au natură pozitivă, iar apoi domeniul lor de aplicare sunt Sisteme pozitive, sau Sisteme care ar putea fi transformate într-un sistem pozitiv echivalent cu o semnificație fizică xor Tirst order delay sistem évolution Timp Figura Evoluția sistemului de întârziere a ordinului setei în modelele constrânse și neconstrânse ■ ■ -■ > -· - JG Concluzii O vedere relaxată a continuilor a sistemelor cu evenimente discrete, a rețelelor Petri continue, a fost luată în considerare împreună cu Diagramele Forrester și sistemele de ecuații diferențiale ordinare liniare (în principal sisteme pozitive) Ele au fost comparate pentru a obține o cunoaștere mai profundă a expresivității relaxării continue a PN-urilor sub semantică infinită a serverului comportamentele pot fi modelate cu PN-uri continue semantice de server infinite Arcurile de control ideale constituie o extensie pentru semantică infinită a serverului PN-uri continue (sunt arcele de control ponderate cu factorul <»), cu ele puterea sa expresivă crește deoarece permit simularea oricărui sistem liniar pozitiv (mărginit sau nu) sau oricare altul care poate fi transformat într-un sistem pozitiv Referencés [ ] M Silva și L Recalde, „Rețele Petri și relaxări integrale: O vedere a rețelelor Petri continue”, IEEE Trans Despre sisteme, om și cibemetică, ( ): - , [ ] E Jiménez, L Recalde și M Silva, „Diagrame Forrester și rețele Petri continue: O vedere comparativă”, În Proc Din cele * IEEE Int Conf On Emerging Technologies and Factory Automation (ETFA ), Nize, pp - , [ ] M Silva și L Recalde, „Unforced Continuous Petri Nets and Positive Systems Sisteme pozitive”, Proc Primul simpozion internațional multidisciplinar privind sistemele pozitive: teorie și aplicații (POSTA ), , pp - , [ ] T Murata, „Rețele Petri: Proprietăți, analize și aplicații”, Proc of the IEEE, ( ): - , [ ] M Silva „Prezentarea rețelelor Petri”, Practice of Petri Nets in Manufacturing, p - Ch & H,, [ ] L Recalde și M Silva, „PN fluidification revisited: Semantics and steady state”, APII-JESA, ( ): - , [ ] H Alla și R David, „Continuous and hybrid Petri nets”, Journal of Circuits, Systems, and Computers, ( ): - , [ ] Jay W Forrester, Industrial Dynamics MIT Press, Mass, Cambridge, [ ] L Farina și S Rinaldi, Sisteme liniare pozitive Teorie și aplicații Matematică pură și aplicată, John Wiley and Sons, New York, Rețelele Petri continue sub semantică infinită a serverului conduc sisteme liniare în parte, prevăzute cu stare și ieșiri nenegative (au o constrângere internă sau implicită de non-negativitate) Arcurile de control ponderate cu factorul k sunt o abreviere în infinité server semantics continuu PN și permit simularea sistemelor liniare pozitive mărginite (eventual sub anumite transformări) Dar comportamentul și puterea expresivă a PN-urilor continue cu semantică infinită a serverului nu sunt limitate la sistemele pozitive liniare mărginite De fapt, orice sistem liniar mărginit (pozitiv sau nu) poate fi mutat către realii pozitivi și, prin urmare, modelat printr-o semantică de server infinită PN continuă În special pur oscilator 